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The influence of a periodic landscape of pinning sites on the vortex dynamics in Pb thin films 
is explored by ac-susceptibility measurements. For different amplitudes h of the ac-drive, the ac- 
susceptibility % — X + *x" exhibits a complex field dependence associated with different dynamic 
regimes. At very low ac-drives where both, multiquanta vortices trapped by the antidots and 
interstitial vortices oscillate inside the pinning potential (intra- valley motion), a small kink in x'(H) 
together with a very low dissipation is observed. At intermediate ac-excitations such that vortices in 
the antidots remain pinned whereas interstices are driven out of the pinning well, a more pronounced 
kink in the screening coinciding with the onset of dissipation (\" (H) 7^ 0) indicates the entrance of 
interstitial vortices. Eventually, at high enough amplitudes all vortices are set in motion and the 
penetration of interstitial vortices appears as a sudden reduction of the shielding power. We show 
that these distinctive features allow us to clearly identify the maximum number of flux quanta per 
hole regardless the vortex dynamic regime. 
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I. INTRODUCTION 

An increasing effort has been devoted in the last years 
to study the interactions of the flux line lattice with a 
periodic pinning landscape This problem represents 
a particular case of a periodic elastic medium interacting 
with an array of obstacles which appears in many other 
physical systemsi Vi 6 

It has been showni that the introduction of a regu- 
lar array of pinning centers leads to strong commensu- 
rability effects of the vortex lattice at matching fields 
H = nHi , where Hi w &o/d 2 , $0 is the superconducting 
flux quantum and d is the period of the pinning array. 
These effects manifest themselves in different ways de- 
pending on the degree of stability of the vortex configu- 
ration and the nature of the vortex pinning. For fields 
H < H s = n s Hi, where n s is the maximum number of 
flux quanta that a pinning center can trap (the satura- 
tion number), every multiquanta ro<I>o-vortex (n < n s ) 
is strongly pinned by an individual defect. A particular 
case arises at H — n s Hi, commonly described as a Mott 
insulator phase with a large compressional elastic modu- 
lus Cn which may give rise to a fixed density of vortices 
over a finite range of HM& 

A different situation occurs for H > n s Hi where fully 
occupied pinning centers are no longer attractive poten- 
tials but interact repulsively with the incoming vortices. 
In contrast to the localized Mott or insulating-like phase, 
in this field range, a metallic-like behavior is observed 
due to the weakly pinned interstitial $o-vortices. The 
peculiar transition at H = n s Hi, from localized to delo- 
calized particles is characterized by a sudden reduction 
of the depinning force density together with a substantial 
increase in the average vortex mobility. For H > n s Hi, 
commensurability effects between the interstitial vortex 
configuration and the pinning array, give rise to local in- 
crements of the critical current J c >I2iiiii£ 



Although these effects have been mainly stud- 
ied by transport and dc-magnetization measurements, 
less explored has been the ac response of these 
systems ^Liiiii^ Ac-susceptibility allows one to readily 
cover a wide range of dynamic regimes from a passive 
linear regime, preserving the initial vortex configuration, 
to the highly disturbing critical state regime, by simply 
sweeping the amplitude h of the ac-drive — 

In the present work we study the field and temperature 
dependence of the vortex pinning in Pb thin films with a 
square array of micro-holes by means of ac-susceptibility 



X = X + 



ix" ■ We found that x{H) exhibits a complex 



behavior associated with a rich variety of vortex dynam- 
ics phases. The possibility to access different regimes 
by tuning the amplitude of the drive field enables us to 
determine the saturation number n s for several temper- 
atures T and dc-fields H, which is difficult to achieve by 
other techniques. 



II. EXPERIMENTAL DETAILS 

The experiments were conducted on two Pb thin films 
with a square antidot array of period d — 1.5 /im, which 
corresponds to a first matching field Hi = 9.2 G. The 
antidots have a square shape with a size b = 0.8 /im. 
The tag, thickness 5 and critical temperature of the used 
samples are, AD100 (5 = 100 nm, T c = 7.22 K) and 
AD65 (6 = 65 nm, T c = 7.21 K), respectively. From 
the T C (H) slope we have estimated, for both samples, a 
superconducting coherence length £(0) ~ 33 ± 1 nm. 

It is important to mention that even though bulk 
Pb is a Type-I superconductor with a Ginzburg-Landau 
parameter— k = 0.45 < l/v2, thin films of this material 
become Type-II superconductors below a certain critit- 
cal thickness S c m& In the particular case of Pb films the 
critical thickness 19,20 is S c ~ 250 nm and therefore, both 
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samples used in the present work exhibit a clear mixed 
state characteristic of the Type-II superconductors. In 
addition, it has been recently shown that this effect be- 
comes more pronounced in patterned samples where the 
effective penetration depth is further increased due to the 
presence of the nano-structuringi^l 

The predefined resist-dot patterns were prepared by 
electron-beam lithography in a polymethyl metacry- 
late/methyl metacrylate (PMMA/MMA) resist bilayer 
covering the Si0 2 substrate. A Ge(20 A)/Pb/Ge(200 A) 
film was then electron-beam evaporated onto this mask 
while keeping the substrate at liquid nitrogen tempera- 
ture. Finally, the resist was removed in a lift-off proce- 
dure in warm acetone. 

The ac-measurements were carried out in a commercial 
Quantum Design-PPMS device with drive field ampli- 
tudes h ranging from 2 mOe to 10 Oe, and the frequency 
/ from 10 Hz to 15 kHz. In this frequency window we 
have found that \ depends only weakly on / and there- 
fore we report results obtained at the same frequency 
/ = 3837 Hz. In all cases, the data were normalized to 
have a total step Ax' = 1, with H =0 at T < 5 K. 
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III. RESULTS AND DISCUSSION 

As it has been demonstrated in previous theoretical^ 
and experimental studies^i^^ rearrangements of the 
vortex lattice at the matching fields nHi have a profound 
effect on the critical current. Accordingly, these changes 
in the pinning properties of the flux line lattice should be 
also reflected in changes of the efficiency to screen out a 
perturbative external ac-field^i In order to corroborate 
this effect we have measured x as a function of H at fixed 
T and h. Figure^a) shows the x'{H) (open circles) and 
x"(H) (filled circles) dependencies for the AD100 sample 
at T = 7.09 K and h = 90 mG^ 

At low H, a smooth decrease of the screening together 
with a very low dissipation x" ~ 0) is observed as H is 
increased. At H = 3iJi, a well defined kink in x'(H) 
coinciding with the onset of dissipation, occurs. This be- 
havior indicates that at H = 3-ffi a substantial change 
in the pinning nature takes place. We argue that at this 
field it is no longer energetically convenient for incoming 
vortices to sit into the antidots and consequently, intersti- 
tial vortices appear. According to this scenario, although 
for H < 3f/i the current induced by the ac-drive is not 
large enough to move out of the potential well vortices 
pinned by the antidots, for H > 3Hi the same force can 
easily set in motion the weaker pinned interstitial vortices 
giving rise to a clear change in the screening properties. 
This behavior is consistent with recent transport mea- 
surements conducted on similar samples £ In brief, the 
sudden reduction of the effective pinning force density at 
H = 3Hi allows for the used ac-amplitude to connect the 
ascending and descending branches of the dc-hysteresis 
loop giving rise to a detectable dissipation^ For lower 
dc-fields, the minor loops built with the ac-cycle do not 



FIG. 1: (a) Field dependence of the real \' (open sym- 
bols) and imaginary x" (solid symbols) components of the 
ac-susceptibility for the AD100 sample at h = 90 mG. (b) 
Field dependence of the critical current extracted from the 
x'(H) curve shown in (a) (see text). The shadowed area in- 
dicates the multiquanta vortex regime. 

bridge the magnetization branches and the dissipation is 
negligible. 

For fields 3Hi < H < 6Hi a series of peaks in the 
screening are seen at each matching field, in agreement 
with the experimentally observed local enhancement of 
the critical current ^SiALi^ These maxima in the screen- 
ing are accompanied by local dips in the dissipation. In- 
terestingly, at 6-ffi the screening undergoes a dramatic 
reduction, which coincides with the maximum dissipa- 
tion. For higher fields, small peaks are still visible at 
7 Hi and 8Hi, before the system enters the normal state. 
Now, these small peaks correspond to local maxima in 
the dissipation. 

As we will demonstrate below, the observed sudden 
reduction at 6Hi can be an artifact owing to the higher 
sensitivity of the screening in that particular range of ac- 
penetration depths (A) and it is not necessarily related 
with a change in the vortex dynamic regime. Indeed, 
for the sake of clarity consider for instance that we can 
approximate our sample shape by a strip with both, H 
and h, applied perpendicular to the plane of the strip. In 
this case, according to the Bean critical state model, the 
components of the ac-susceptibility are given by^ 

x' = -^, (1) 
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// _ Tanh(s) . Tanh(s/2) / n \ 

X i 1 JJ2 > \ z ) 

where s w A/S oc J^ 1 . The derivative of the eq.QJ and 
therefore the sensitivity of x' to small changes in A, max- 
imizes at s = 1, i.e. when the Bean penetration depth 
coincides with the sample thickness. On the other hand, 
the situation A ~ 6 is achieved at the maximum of the 
dissipation x"r^ thus naturally accounting for the sys- 
tematic coincidence of the field position of the maximum 
dissipation and the jump in the screening. Furthermore, 
using the strip geometry approximation we can estimate 
from x{H) the field dependence of the critical current J c . 
The result of this analysis is shown in Figure^b). In this 
Figure we observe that the height of the jump at 6Hi has 
now no particularities and is similar to that observed at 
4H"i. In contrast to that, the distinct feature signaling 
the entrance of interstitial vortices at 3Hi is still clearly 
seen. 

It is important to note that the kink in x'(-ff) at 
H = 2>H\ allows us to determine the saturation number 
n s (T, b, £), which is an intrinsic property of the pinning 
array and therefore it is independent of h. In contrast to 
that, the big jump associated with the peak in the dissi- 
pation, should occur at a field H such that J C (H) w h/6 
and thus should depend on the strength of the ac-drive 
h. This is in fact confirmed by the x(H) measurements 
performed on the same sample at the same temperature, 
with a higher excitation h = 500 mG (see Figure EI a ))- 
In this case, we observe that the largest reduction of the 
screening is shifted down to H s = 3-Hi following the po- 
sition of the maximum dissipation. 

Performing the same procedure described above, we es- 
timate the critical current J C (H) obtained from the x'{H) 
curve at h = 500 mG (see filled circles in Figure H£b)). 
Now we can see that after this transformation the sharp 
change at H s — 3Hi persists and is not an artifact due to 
the amplification in the screening sensitivity. Since, ac- 
cording to the Bean model J c is independent of the used 
ac-excitation, the critical currents obtained for h = 90 
mG and h = 500 mG should coincide in the field range 
where this approximation holds. This behavior becomes 
apparent in Fig|2Jb) where we superimpose the J C (H) 
curve for h — 90 mG, already shown in Fig^b) (open 
circles) . The remarkable coincidence of the data for fields 
H > H s indicates that in this particular field range the 
simple Bean model is a very robust approximation and 
correctly accounts for the observed behavior. On the 
other hand, the difference observed at H < H s is due 
to the fact that for h = 90 mG the used Bean critical 
state is no longer valid and therefore the J C (H) curve so 
determined is ill-defined. 

Another feature shown in FigJ2J is that for H < H s , 
where the antidots act as attractive centers, commensu- 
rability effects manifest themselves as small reductions 
in the screening as H increases. In this field range, the 
vortex-antidot interaction dominates over vortex-vortex 
interaction and thus, the decrease in the effective pinning 
potential strength at nli\ will result in a reduction in the 



shielding power. For H > H s , filled antidots become re- 
pulsive centers and an opposite behavior is observed^ At 
these high fields, vortex- vortex interactions are more rel- 
evant and consequently, collective effects give rise to local 
increments in the screening properties at every commen- 
surability field. We have verified that these features sig- 
naling the transition from attractive to repulsive vortex- 
pin site interactions which allow one to identify the satu- 
ration number n s , were present in six other samples with 
periodic pinning array. 
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FIG. 2: (a) Field dependence of the real \' (open sym- 
bols) and imaginary \" (solid symbols) components of the 
ac-susceptibility for the AD100 sample for h = 500 mG, at 
t = 0.980. (b) Critical current as a function of field, obtained 
from the x'{H) curves for h — 90 mG and h = 500 mG, at 
t = 0.980 (see text). The shadowed area indicates the multi- 
quanta vortex regime. 

Let us now analyze the experimentally obtained n s = 
3 value. According to the theoretical estimation of 
Mkrtchyan and Schmidt (MS) r& the number of trapped 
vortices by a hole is n s ss 6/4£(T) which in our partic- 
ular case gives n s ~ 1. Although this value is smaller 
than the experimentally determined above, it is impor- 
tant to note that this model underestimates the real n s 
value since it considers a single cylindrical cavity with 
radius b <C A. Clearly, this hypothesis is not fulfilled 
in our system where b > 2X(T). A more detailed anal- 
isys of the validity of this model could be obtained from 
the temperature dependence n s (T). However, matching 
effects reported here are restricted to a very narrow tem- 
perature ranged close to T c where n s remains almost 
constant and no evident temperature dependence is ob- 
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served. Additionally, it is also expected that n s increases 
as the applied de-field is increasedS^Si^i An extension of 
the original work of Mkrtchyan and Schmidt to arbitrary 
large cavity radius has been recently done by Nordborg 
and Vinokur— using the London approximation. Ad- 
diotonally, by simple energetic considerations Buzdin 31 
showed that in a triangular vortex lattice, two-quanta 
vortex becomes energetically favorable for temperatures 
such that b 3 < £(T)A(T) 2 , a condition that, m our sam- 
ple, is satisfied for t < 0.995. The experimental evidence 
on this point also suggests that a higher filling than that 
predicted by MS has to be considered. Indeed, Bitter 
decoration experiments performed by Bezryadin et aljS 
on a Nb film with a triangular array (period d = 6.1/j,m) 
of blind holes (size b — 0.8/im), showed that n s = 3 very 
close to the critical temperature for H — 6.37 G. More 
recently, Grigorenko et ali^4 using scanning Hall probe 
microscopy in a Pb film with a square array of antidots 
of period 1.5/zm and hole size b = 0.66/im showed that 
n s = 2 at t = 0.77. Although this number is some- 
how smaller than our estimate, it should be noted that 
in those experiments, the measurements are performed 
in a field-cooled metastable state which corresponds to 
a frozen vortex configuration at a higher temperature. 
On top of that, since the n s (b) is a stepwise increasing 
function^ even though the holes used in their work are 
slightly smaller than ours, this small difference can eas- 
ily lead to an extra vortex trapped per hole. Thus, in 
general, our determination of the saturation number by 
means of ac-susceptibility is consistent with the previ- 
ously reported results using alternative techniques. 

So far we have shown that at low h amplitudes, the cur- 
rent induced by the ac-drive is not large enough to move 
vortices trapped by antidots out of the pinning poten- 
tial but can set in motion the weakly pinned interstitial 
vortices. According to this picture, it should be possible 
to reach a situation where interstitial vortices perform 
intra- valley oscillations as well, by further decreasing the 
ac-drive. To check this possibility we have measured the 
ac response of the sample AD65 at very low amplitudes 
(h = 3 mG and h = 6.4 mG), as shown in the main panel 
of Figure The dotted line indicates the upper limit 
of the field range (H < AH\) where the detected dissi- 
pation lies below our experimental resolution (x" ~ 0) 
and therefore where vortices oscillate inside the pinning 
potential. Within this regime the ac-response is inde- 
pendent of h, confirming that the system is in the linear 
regime. We also note that a clear kink at H = 3Hi 
signals the transition from a intra-valley motion of mul- 
tiquanta vortices trapped by the antidots, to an intra- 
valley motion of interstitial vortices inside the potential 
produced by the strongly pinned neighbors and vortices 
trapped by the holes. For fields H > AH\ a finite dis- 
sipation together with an amplitude dependent response 
indicate the onset of inter-valley motion of interstitial 
vortices, whereas vortices in the antidots remain pinned. 

We turn now to the analysis of the field dependence 
of dissipation x" ■ As we noticed previously, for fields 
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FIG. 3: Main panel: x(H/Hi) for the AD65 sample for h = 3 
mG and h = 6.4 mG, at t = 0.988. The dotted line indi- 
cates the onset of dissipation. Inset: x{H/Hi) for the AD100 
sample at h — 500 mG and t — 0.97. 



smaller than that where the maximum dissipation takes 
place, the local maxima in the screening, and therefore in 
the critical current, are accompanied by local dips in the 
dissipation. In contrast, for higher fields, peaks in \x'\ 
give rise to local increments of x" ■ This effect is more 
clearly seen in the inset of FigEI where we show x{H) 
for the AD100 sample at t = 0.97 and h = 500 mG. 
The observed behavior can be easily understood within 
a critical state scenario. Indeed, in the simplest case 
of a Bean critical state model for a thin sample with 
applied field H perpendicular to the plane of the film (see 
eq.Q and eq.Q), the maximum dissipation occurs when 
h reaches the full penetration field H d (T,H) = J e <5i 25 i 26 
According to this model, the asymptotic limits of the real 
and imaginary ac-susceptibility components are related 
by, 



(i-lxU 



h < J C S 



(3) 



X 



X 



h > J c 5 



From these equations we can see that for low dc-fields 
(h -C J c 5) an enhancement of the screening leads to a 
suppression of the dissipation, whereas an opposite be- 
havior is predicted for high fields (h 3> J c 5), in agreement 
with our experimental observation. It is interesting to 
note that in this case the commensurability effects have 
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the benefit of acting as markers in J(H) which eventually 
help to study and identify the different dynamic regimes. 

Finally we would like to note that no clear indications 
of the depinning of vortices trapped by the antidots have 
been observed. This depinning process of multiquanta 
vortices at H t represents a peculiar transition from a 
state where supercurrents flow around an antidot (for 
H < H t ) with the n$o-vortex to a n depinned <&o-vortices 
state (for H > H t ). This metamorphoses of single to mul- 
tiple entities is a very interesting issue that has not been 
addressed so far and clearly deserves further theoretical 
and experimental investigations. 



a reliable determination of the saturation number n s of 
the pinning structure which is hard to find by other indi- 
rect methods. Second, at a higher H the enhanced sen- 
sitivity of the ac-screening, when the penetration depth 
A <~ t, accents the drop on the critical current leading 
to a dramatic reduction of x'(H). In general, this anal- 
ysis represents a powerful method to study the strength 
of different pinning centers arranged periodically and in 
particular, provides an alternative way for determining 
the maximum number of flux lines trapped per pinning 
site as a function of temperature. 



IV. CONCLUSION 

In summary we have shown that the field dependence 
of the ac response of samples with a periodic pinning ar- 
ray exhibits a rich behavior associated with different dy- 
namic regimes. In particular, we show that for low drive 
fields two distinctive features emerge. First, at low dc- 
fields H, a kink in the screening x'{H) indicates the en- 
trance of interstitial vortices in the sample. This provides 
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